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■ 1. Introduction 



In this work we want to explore the relationship between certain eigenvalue condition 
for the symbols of first order partial differential operators describing evolution processes 



and the linear and nonlinear stability of their stationary solutions. 



Consider the initial value problem for the following general first order quasi-linear 
system of equations 

d 

v t = P(v, x, t, V)i> = > A u (v, x, t)— — v + B(v, x, t)v, 

~ ox v 

v(x,0) = f(x). 

Here v is a (column) vector valued function of the real space variables (xi, . . . , x s ) and 
time t with components v±, . . . , v n . A v and B are nxn matrices and f(x) is a vector valued 
function of the space variables. 



'* Work supported by the Office of Naval Research n00014 90 j 1382 
* Work partially supported by CONICET. 
t Researcher of CONICET. 



1 



We are interested in solutions which are 27r-periodic in all space variables. There is 
no difficulty to extend the results to the Cauchy problem on the whole x-space. Instead of 
Fourier series we would use Fourier integrals. 

We shall restrict our considerations to the case 

u t = ^2(A Qu + eA lv (x,t, u, a))- — u+ (B + eB 1 (x,t,u,e))u. (1.1) 

u=0 ° Xv 

Here Aq u , Bq are constant matrices and £ is a small parameter. This is, for instance, the 
case when the stationary solution is constant and and we consider the solution close to the 
steady state. 

Assumption 1.1. For every p = 0, 1, 2, . . . and any c > 0, there is a constant K p such 
that the maximum norm of the p th derivatives of A\ v , B\ with respect to x,t,e and u are 
bounded by K p , provided \u\qq < c. For f(x), the corresponding estimates hold. 

Definition 1.1. The system (1.1) is said to satisfy the stability eigenvalue condition 

if there is a constant 5 > such that, for all real uj, the eigenvalues A of the symbol 

s 

P (iuj) + B :=iJ2 A o^ u + Bo (1.2) 

satisfy 

ReA<-5. (1.3) 
We have to define stability for system (1.1). 

Definition 1.2. The system (1.1) is stable if, for any /, there exists an eq such that, for 
< e < Eq, the solutions of (1.1) converge to zero for t — > oo; and there exists an integer 
Po such that eo depends only on the constants K p with p < p .l 

In this work we shall look at sufficient conditions under which the stability eigenvalue 
condition implies stability. 

Consider first the constant coefficient case, i.e., set e = in the above system. In 
Section 2 we shall prove that it is possible to find a positive definite selfadjoint operator 
Hq such that all solutions of the system satisfy 

^(u,H u) < -8(u,H u), 

provided that the problem is well posed in the Li sense and the eigenvalue condition is 
satisfied. In this case the system of equations is a contraction in a new norm. 

In Section 3 we consider linear systems with variable coefficients, i.e., the A lv depend 
on x and t but not on u. The construction of H proceeds via the theory of pseudo-differential 
operators, i.e., we construct the symbol H(x,t,u>) and define the operator H by 

Hu = Y,e l{LJ ' x) H(x,t,uj)u(uj) for all u = ^ e l ^^u(u). 

1 We have not specified the norm under which that convergence takes place, but we 
shall be using uniform pointwise convergence. 
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H depends on the symbols 

s s 

Po(iu) =: ' s ^2 / Aq v ioo v , Pi(x,t,iu) =: ^2A lu (x,t)iu> l/ . 

We need that if is a smooth function of all variables. This is only the case if P , Pi satisfy 
extra restrictions. For the linear and the nonlinear case, we make one of the following 
assumptions. 

Assumption 1.2. The stability eigenvalue condition is satisfied and the multiplicities of 
the eigenvalues of Po(iu>) + sPi(x,t,u,iui) do not depend on x, t, u, u, e. Also, for every 
t, u, u>, s, there is a complete system of eigenvectors. 

Assumption 1.3. The stability eigenvalue condition is satisfied and the matrices Aq u , Bq 
and Aiv, v = 1, . . . , s, are Hermitian. 

Under any of these conditions we can again construct an if-norm and prove that the 
problem becomes a contraction. 

In the last section we consider the nonlinear equations and the main result of this 
paper is 

Main theorem. Suppose that Assumption 1.1 and Assumption 1.2 or 1.3 hold. Then, 
for sufficiently small e, the problem is a contraction in a suitable H-norm and the system 
(1.1) is thus stable. 

In the Appendix we relax the eigenvalue condition somewhat. 

To prove stability for time dependent partial differential equations via changing the 
norm has been done before. For example, in [1] the method was applied to mixed symmetric 
hyperbolic-parabolic equations which included the Navier-Stokes equations. In that case 
H was explicitly constructed and not related to an eigenvalue condition. If we make 
Assumption 1.3, then our H is similar to the H in [1]. 

2. Systems with constant coefficients 

In this section we consider the system 

* = E^|£ + ^= : + (2i) 

y(x,0) = f(x), 

with constant coefficients. We are interested in solutions which are 27r-periodic in all space 
variables. We assume that the problem is well posed in the Li sense, i.e., for every T there 
exists a constant K(T) such that the solutions of (2.1) satisfy the estimate 

\\y(;t)\\<K(T)\\y(;0)l < t < T. (2.2) 

Here 

(u,v) = ' ■ '/ (u,v)dxi ■ ■ ■ dx s , \\u\\ 2 = (u, u), 
Jo Jo 
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denote the usual L 2 scalar product and norm. 

One can characterize well posed problems algebraically. Using the Kreiss matrix 
theorem (see [2], Sec. 2. 3), one can prove 

Theorem 2.1. The problem (2.1) is well posed in the L 2 sense if and only if it is strongly 
hyperbolic, i.e., the eigenvalues of the symbol 

s 

Po(iuj) = i Aq v oo v , ujj real, 

are purely imaginary and, for every fixed u>' = oj/\oj\, there exists a complete set of eigen- 
vectors ti, . . . , t n which is uniformly independent, i.e., there is a constant K such that 

\T~ l \ + \T\<K, T =(*!,..., t n ). 



We can expand the solution of (2.1) into a Fourier series 

y(M) = Xy< w '*>y(u/,f). (2.3) 

The Fourier coefficients are the solution of the Fourier transformed system (2.1) 

s 

y t =(iJ2 A „uj u + S ) y =: (PoM + So) V- (2.4) 

We assume that the eigenvalue condition (1.2), (1.3) is satisfied. Then we can find, for 
every fixed u>, a positive definite Hermitian matrix H, a Lyapunov function, such that 

2ReH(cj)(P (iuj) + Bo) =: H{u)(P {iu) + B ) + (P * (iu) + B *) H (u) 

(2.5) 

< -8H{u). 

Therefore, 

— t), H(u)y(u, t)) = 2Re(y(iv, t),H(u) (P (iu) + B )y(u, t)) 

< -5{y{u,t),H{u)y{u,t)). 

Thus, for every fixed ui, the transformed system (2.4) is a contraction in the (a;)-norm. 
Using the Kreiss matrix theorem, one can prove (see [2, Sec. 2. 3]) 

Theorem 2.2. Assume that the problem (2.1) is well posed in the L 2 sense and that the 
eigenvalue condition (1.2), (1.3) is satisfied. Then we construct the matrices H(ui) such 
that they satisfy the uniform inequalities 

Kl x I < H{uj) < KJ. (2.6) 
4 



Here K 4 does not depend on u>. 

We can use H(u) to define an operator H by 

Hu = J2 H(u)u(u)e l{uj ' x l (2.7) 

UJ 

It has the following properties 

(1) H is selfadjoint and 

K^Wuf < (u,Hu) < K 4 \\u\\ 2 . 

(2) 2Re H(P + B ) =: H(P + B ) + (P * + B^)H < -SH. 
These properties follow from Parseval's relation 

uj 

= Y j (H{uj)v{u)Mu)) = {Hv,u). 

uj 

Also, 

K^\\uf = Kl 1 H^)\ 2 < 5>(u/)>£(u/)«(u/)> 
= (u,Hu) < K 4 \\u\\ 2 

and 

2(u,ReH(P + B )u) = 2^( 

UJ 

< -5^{u(uj),H(uj)u(uj)) = -5(u,Hu). 

UJ 

Thus, we can use H to define a new scalar product by 

(v,u) H = (v,Hu), \\u\\ 2 H = (u,u) H , 
which is equivalent with the L2-norm. The second property gives us 

Theorem 2.3. If the conditions of Theorem 2.2 are satisfied, then the problem (2.1) is a 
contraction in the H-norm. 

Proof. 

-(y, Hy) = 2Re(y, H(P + B )y) < -8(y, Hy). 
This proves the theorem. 

3. Linear systems with variable coefficients 
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In this section we want to generalize Theorem 2.2 to linear systems 



s 



v t = J2( A °» + £A i»( x i *)) jr~ + ( B o + eB i) v 

»=o ° Xv (3.1) 

and show that it is a contraction in a suitable if-norm. We shall construct the if-norm 
with help of a pseudo-differential operator 

H{t) = H ^S^eH 1 {t) (3.2) 

with the following properties. 

(1) Hq, S, Hi(t) are bounded selfadjoint operators. Hq and S do not depend on t. dH\/dt 
exists and is also a bounded operator. Thus, there is a constant K such that 

||ffo|| + ||5|| + ||#i(*)|| + \\^\\<K. 

(2) Hq + S is positive definite with K such that 

\\H Q + S\\ + \\(H + sy^i <K. 

(3) 2Re H P =: H P + P * #o = 0. 

(4) 2Re(tf + S)(P + B ) = 2Re(SP + H B ) < -S(H + S). 

(5) S is a smoothing operator with 

||SPi|| < K. 

(6) \\Re(H + eH l (t))(Po + eP 1 )\\=e\\Re(H P 1 + H 1 P + eH 1 P l )\\ <eK. 
We can prove 

Theorem 3.1. Assume that there is an operator H of the form (3.2) with the proper- 
ties (l)-(6). For sufficiently small e the scalar product (u, Hv) defines a norm which is 
equivalent with the L2-norm and the system (3.1) is a contraction in the H-norm. 

Proof. That (u, Hv) defines a norm which is equivalent with the L2-norm follows from 
properties (1) and (2). Also, 

— (u, Hu) = s(u, H lt u) + 2Re(u, (H + S + sH 1 )(P + B Q + e(P x + Bi))u) 

= s(u, H lt u) + 2Re(w, (H + S){P + B )u) 
+ 2Re(w, (H + eH 1 )(P + eP^u) 

+ 2eRe(u, {H B X + 5(Pi + B x ) + H^Bq + eSi))u) 

< -(5 + 0(e))(u,Hu). 



This proves the theorem. 

We construct the symbol of the pseudo-differential operator (3.2) in the following way. 
Consider all systems with constant coefficients which we obtain by freezing the coefficients 
of (3.1) at every point x = xo, t = to. We assume that the initial value problem for all 
these systems is well posed in the L>2 sense and, therefore, we can construct the matrices 
H(x, t, ui) for every fixed x, t. Now we think of H(x, t, u) as a symbol of a pseudo-differential 
operator where x, t are independent variables. Formally, we define the operator H by 

LO 

This definition makes sense only if H satisfies the usual properties of symbols for pseudo- 
differential operators. Also, we need the algebra for such operators to prove that (3.1) 
becomes a contraction. We want to prove 

Theorem 3.2. Assume that the following conditions hold. 

a) There exists a positive definite Hermitian matrix Hq(uj') which is a smooth function of 
lo' = u>/\u>\ such that 

2Re H (uj')Po(ioj) =: H (aj')P (iuj) + Pq(iuj)Ho(oj') = 0. (3.3) 

b) For sufficiently large \u\, there is a Hermitian matrix S = S(oj', 1/\uj\) which is a smooth 
function of u' and such that 

2Re(tf (^) + 7^S(u', {\u\P (iu') + B ) < -S(H (u') + t^S(u', (3.4) 

r 4 ' 1^1 

c) There exists a Hermitian matrix Hi(x, t, u>') which is a smooth function of x, t, u>' such 
that 

2Eb(h (u') +eH 1 (x,t,u',ej^ (P (iu') + eP 1 (x,t,iu'j S ) =0. (3.5) 

Then we can construct the pseudo-differential operator (3.2) which has the properties 
(l)-(6). Also, there exists an integer po such that the constant K depends only on the 
first po derivatives of the symbols and of the coefficients of (3.1). Thus, the problem (3.1) 
is a contraction in the H-norm. 

Proof. We construct the symbols for the pseudo-differential operators 

H u = J2 et(uJ,x) Ho(oo)u(uj), 

LO 

Su = Y j e l ^ x) S{u)u(u>) (3 . 6) 

LO 

H lU = J2 ei{uJ,x) Hi{x,t,u)u{iv). 

LO 

Ho(u>), S(u) do not depend on x,t. 
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Let C > be a constant. Consider the symbol (1.2) for \uj\ < C. The inequality (1.3) 
implies (see Lemma 3.2.9 in [2]) that there is a positive definite Hermitian matrix S^'{u) 
which is a smooth function of u> such that 

2ReS {l \u)[P {iu) + B ) <-8S {1) (oj), \w\<C + l. 
Let </?(|cu|) G C°° be a monotone cut-off function with 

v(M) = 



1 for|u;|>C+l 
for U < C 



We define 



H q {uj) = ^{\uj\)H q {u'), 
S(u) = V^SV, 1/M) + (1 - V (|a;|))5W(a;). 

It follows from (3.3) and (3.4) that, for sufficiently large C, the operators H and S 1 have 
the properties (l)-(5). The symbol 



( p(\(j\)(H ((j') + eH 1 ((j , t x t t)) 



defines a pseudo-differential operator Ho + eH\\ and the algebra of such operators shows 
that 

H + eH x =H + |(#n + H* lx ) (3.7) 

has the desired properties (1) and (6) and K can be estimated as required. This proves 
the theorem. 

We shall now give algebraic conditions such that the conditions of Theorem 3.2 are 
satisfied. 

Theorem 3.3. Assume that Assumption 1.2 holds. Then we can construct the symbols 
of Theorem 3.2 whose derivatives can be estimated in terms of the derivatives of the 
coefficients of (3.1). Therefore, for sufficiently small e, the system (3.1) is a contraction. 

Proof. We consider the symbol Pq(iuj) + B = \uj\Pq{iuj') + B in a neighborhood of a 
point u'q. Let Ai, . . . , A r denote the distinct eigenvalues of P (iuj'). It is well known (see, 
for example [3]) that, because of the constancy of the multiplicity of the eigenvalues of 
Po(iu'), there exists a smooth nonsingular transformation Tq(u') such that 

f - 1 (u')P (tuj')f (u')= ^ •.. |. (3.8a) 

All eigenvalues of Aj are equal to Xj and, since there is a complete set of eigenvectors, 



is diagonal. 

T is not unique. We can replace it by 



7o=T 







(3.86) 



Here the T j denote arbitrary nonsingular submatrices. We shall choose them as constant 
matrices later. (3.8a) gives 

T -^|Po(*cy) + Bo)To=: 



'Ai 



\ //i n 
\B r i 



■ ■ B lr 



and (3.8b) gives 

TQ\\u\P Q {iu') + B )T Q 



Ai 



/ Tqx BuTqi B 



\U!\ 



+ 



Ar 



12 



B 



lr 



\ 



(3.9) 



\ B r \ B r 2 ■ ■ ■ T 0r , 1 5 rr Tor / 



For large \oj\, we can consider the second matrix in (3.9) as a small perturbation of the 
first. Therefore, (again, see [2]) there is a smooth transformation T\(uj' ', l/|w|) such that 

(I + ^T 1 )- 1 T - 1 (| W |Po(^^ , ) + B )T (I + -LrO 

\uj\ \UJ\ 



'Ax 



= M 



+ 



A T 



\ 



V 



o 



Tq t B rr To r -\- ,,B rr J 



By assumption, the eigenvalues of Bjj have negative real parts. Therefore, we can choose 
T j such that 

38 

2Re(T~ 1 B jj T 0j ) < i", \J - u' Q \ sufficiently small. 



(Again, see Lemma 3.2.9 in [2].) Thus, 
and, for sufficiently large \u\, 

^ + R 5 =(( (7 + R ri)ro )"T( (7+ R ri)r °)" 1 

satisfies (3.3) and (3.4). By the usual partition of unity argument, we can construct H 
and S for all u' and conditions (a) and (b) in Theorem 3.2 hold. 
We now consider the matrix (symbol) 

P (iu') + ePi{x,t,iu'). (3.10) 

As the eigenvalues of (3.10) are purely imaginary and their multiplicity does not change, 
we can find a smooth transformation T 2 (x, t, ui', e) such that 

(/ + eT^T^iPo + eA)T (/ + eT 2 ) 




Here Aj = Xjl and T 2 is a smooth function of all variables. The matrix 

H + eH x = (T~ 1 )*T~ 1 , T = T (I + eT 2 ), 

has the property (3.5) and condition (c) in Theorem 3.2 hold. Therefore, Theorem 3.3 
follows from Theorem 3.2. 

We consider now the symmetric systems (3.1), i.e., those satisfying Assumption 1.3. 
In this case the stability eigenvalue condition, for u = 0, implies that 

ReS < -51, (3.11) 

and therefore 

Re(it, (P + B )u) <-5(u,u). 

Thus, we can show that (3.1) is a contraction in the usual L2-norm (H = I). In the 
Appendix we shall relax the eigenvalue condition to some cases where (3.11) does not 
hold. Therefore we give here a proof which does not depend on (3.11). 

Theorem 3.4. Assume that the coefficients Aqj, Aij, j = 1,2, ...,s, and B but not 
necessarily E>i are Hermitian matrices. Assume also that the eigenvalue condition (1.3) 
holds. Then, the results of Theorem 3.3 are valid. 
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Before we give a proof of the last theorem, we will prove 

Theorem 3.5. Assume that, for sufficiently large \uj\, there is a Hermitian matrix H{uj) = 
I + j^S where S = S(uj', 1/M) is a smooth function of to' and 1/\oj\ such that 

2ReH(uj){\u\P (iu') + B ) < -SH(u). (3.11) 
Then, for sufficiently small e, the system (3.1) is a contraction. 

Proof. The proof proceeds similarly to the proof of Theorem 3.2. It is much simpler, 
because in this case we construct a time independent pseudo-differential operator of the 
form 

H = I + S 
which has the properties of Theorem 3.1. 

Proof of Theorem 3.4- Consider the symbol |a;|Po(i^') + Bq for large \u\. Let u>' = uj' be 
fixed. Since the coefficients Aoj are Hermitian, there is a unitary transformation such that 

U*(u' )(\u\P (iu' ) + B ^U(u' ) 




(3.12) 



Here 

A, = Xjl 

represent the different eigenvalues according to their multiplicity. Since Bjj are also Her- 
mitian, we can assume that they are diagonal. Otherwise, we apply a block-diagonal 
unitary transformation to (3.12). For large \ou\, we consider the S-matrix in (3.12) a small 
perturbation of i|u;|A. Therefore, we can construct a transformation 1+ j^T(uj' ) such that 

(/ + pTK)) U*(u' ) (|w|PoK) + B )u(u>' ) (/ + pT^)) 
Mi \ /B n \ 

\ A r / \ B rr ) 

=: i\u\A + B + 7—rB. 

\u>\ 

The eigenvalue condition guarantees that Bjj < —51 for all j and, for sufficiently large |o;|, 

1 ~ 3 
2Re(z|w|A + B+ —B) < --5I. 

\u\ 2 
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We shall now show that there is a neighborhood of uj' where the matrix H(uj) of (3.11) is 
given by 

H{u) = U{u' Q ){l + j— |T*(o;o))~ 1 (I + ^T(^))-V*K) =: / + ±-S{ u ' , -U 

\U)\ \U}\ \U)\ \U)\ 

We have 

2Re H(u){\uj\P (iu') + B ) 

= 2Re H(uj')(\uj\Po(iuj' ) + B ) + \co\2Re H(oj)P (i(uj' - u' )) 

< --SH(u) + \u\ ■ -^-2Re S(iv' , -^)P (i(u' - u' )) 
|cu| \u)\ 

3 



< ^— -5 + const. \oj' — uj' \ ^JH (uj). 



Thus, for sufficiently small \u' — u>' \, the inequality (3.11) holds. With help of the usual 
partition of unity argument (see again Lemma 3.2.9 of [2]), we can construct H(u) for all 
to' and the theorem follows from Theorem 3.5. 

4. Nonlinear systems. 

In this section we consider the nonlinear system (1.1). We start with the case that 
Aq u , Ai„, v = 1, . . . , s; are Hermitian matrices and 

ReB <-5. (4.1) 

Our arguments follow closely the arguments in [2, Chapter 5,6] and we assume that the 
readers are familiar with them. 

We shall derive a priori estimates and shall use the following notations: j = (ji, . . . , j s ), j u 
natural numbers, denotes a multi-index, \j\ = J2j u , -D J = d^/dx^ 1 ■ ■ ■ d js /dx J s s denote 
the space derivatives and 

^=£ii^ii 2 

iii<p 



denotes the derivative norm of order p. 

To begin with, we assume that e = and derive estimates for 



a i=( P °^ +B °>> (4.2) 



u(x,0) = l(x). 

Differentiating (4.2) gives us 



(Diu) t = P (|-)D% + B Diu. 
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Therefore, by (4.1), 

— \\D j u\\ 2 = 2Re( D j u,P (^-)D j u) + 2Re( D j u, B D J u) 
at v ox 

= 2Re(D J u,B D j u) < -2S\\D j u\\ 2 . 
Adding these inequalities for all j with \j\ < p we obtain, for any p, 

— 

di 

i.e., 



< -2S\\u\\l, (4.3) 



IK,0li;<e-^||«(.,0)||J. 

Now we consider the nonlinear system (1.1). We derive an estimate for p > s + 2. Local 
existence causes no difficulty, it has been known for a long time. There exists an interval 
< t < T, T > 0, where the solution exists and 

IK,t)||J<2||«(.,0)||J. (4.4) 

There are two possibilities: 

Either T = oo or T < oo and ||u(.,T)||J = 2\\u(-, 0)||J. (4.5) 

We shall now prove that T = oo for sufficiently small £o an d that the initial value problem 
is a contraction (see [2, Section 6.4.1]) 
We differentiate (1.1) and obtain 

(Di U ) t = (Po(fa) +B )(D3 U ) + eP 1 {x,t,u,— ){D^u) +eRj, 

where Rj denote lower order terms. Therefore, 

d 



\D^uf t = 2Re(l)%, (Po(fa ) + B )iy 



u 



+ 2Re e(p j u, P^x, t, u, — )D 3 i^j + 2eRe(D J u, Rj). 
By (4.1), 

2Re(D j u, (Po(q^) + B Q ^D J v)j < -25\\D J uf. (4.7) 
Integration by parts gives us 



(4.6) 



< const. K 1 (l + ± ^ )\\D^uf 



(4.8) 



< M i\\ u \\l- 
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The Mj are polynomials in \\u\\ p of degree \j\ whose coefficients depend only on the con- 
stants K , . . ■ ,Kj of Assumption 1.1. Using Sobolev inequalities we find bounds 

11(^,^)11 <^H2. 

Adding all these inequalities gives us 

d 

— \\u\\l<-28\\u\\l + eM\\u\\ 2 p , M = maxM J . (4.9) 
Thus, for all e with < e < s$ with s sufficiently small, we have 

J>ll*<-*ll< 

Therefore, T = oo and the initial value problem is a contraction. 

We now consider the general case. We assume that the assumptions of Theorem 3.3 
or 3.4 are satisfied. Again, we begin with the case that e = 0. Then there is a pseudo- 
differential operator 

H = H + S 

which defines a norm that is equivalent with the L2-norm such that 
Thus, 

IK-,t)||^<e- 5t H-,o)||^ 

and, therefore, also 

\H-m, P <e- st \\<M\% p , \wwi P =j:\\Diu\\%- 

\j\<p 

Now we consider the general system (1.1). We proceed in the same way as in the 
previous case and derive estimates for p > s + 2. The only difference is that we derive the 
estimates in the if -norm. 

Local existence is again no difficulty. Thus, there is an interval 0<t<T, T > 
where 

|2 ^ oIL,/ nMl2 



K,t)||^<2|K,o)||^. 



For T, the alternative (4.5) holds. In this interval we can estimate the solution and its 
derivatives up to order p — [s/2] — 1 in the maximum norm in terms of \\u(-, 0) ||^ . Thus, 

we can think of the system (1.1) as a linear system and construct the pseudo-differential 
operator (3.2) and estimate the solution and its derivatives in the i7-norm which differs 
from the i/"-norm only by terms of order s. The symbol depends on the solution but, 
by Theorem 3.2, if p is sufficiently large, then the constant K in Theorem 3.1 can also be 
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estimated in terms of \\u(-, 0) ||^ . The rest of the proof proceeds as before. We differentiate 
(1.1) with respect to the space derivatives and obtain in the if- norm 

— (D J u, HD J u) = e(u, H u u) + (d 3 u,H({P g + B + e{P 1 + B 1 ))D J U S j 
+ 2eRe(D J u,HR J ). 

Using Theorem 3.1, we obtain the inequality (4.9) but now in the if-norm. Thus, we have 
proved the Main Theorem of Section 1. 

Appendix 

We want to relax the stability eigenvalue condition for the cases when some of the 
eigenvalues of So have zero real part. We do this only for the 27r-periodic case. 

Definition A.l. The system (1.1) is said to satisfy the relaxed stability eigenvalue 
condition if the following conditions hold. 

1 ) There is a constant 5 > such that, the eigenvalues A of the symbol P(iuS) + B satisfy 

Re\<-5 (A.l) 

for all u> = (u>i, . . . , uj s ) ^ 0, ujj integer. 

2) The eigenvalues A(0) of B satisfy 

Either Re\<-5 or A = 0. (A.2) 

Also, if the multiplicity of the zero eigenvalue is r, then there are r linearly independent 
eigenvectors connected with A = 0. 

3) The nullsapce of B\ contains the nullspace of Bq. 

We can find a nonsingular transformation S such that 

S^BqS = ^ S 01 jj^, B Q1 nonsingular. (A3) 

If B is symmetric, we can choose S to be unitary. Therefore, we can assume that B 
already has the form (A. 3). Then, by the third part of the assumption, Bi has the form 



Let 



Here the partition of u(0,t) corresponds to that of B ,Bi. Denote by Q the projection 





<MM)=( fl // (M) 



15 



Using the notation 



Qu(x,t) =: w (0) (^), {I -Q)u(x,t) =: v(x,t), 
we can write the system (1.1) as 



i 0) = Q(P + S )(« (0) + v) + eQ(P 1 + B 1 )(u^ + v) 
--eQ(P 1 + B 1 )v. 



v t = (I- Q)(P + B )(« (0) + v) + e(I - Q)(P 1 + B^u™ + v) 
= (P + B )v + e(I - Q)(P 1 + B x )v. 



(A.5) 



(A.6) 



As before, we need only to consider linear systems. Then (A. 6) decouples completely from 
(A.5). It is a system on the subspace (/ — Q)L,2- Our results tell us that, for sufficiently 
small s, it is a contraction and v converges exponentially to zero. Since 

u i0 \t)=e [ Q(Pi + PiMz,0^ W 0) (0), 
Jo 

it follows that also v,(°'(t) converges for t — > oo. 

We summarize the results of the appendix in the following theorem. 

Theorem A.l. Suppose that assumption 1.1 and assumption 1.2 or assumption 1.3 hold 
but with the stability eigenvalue condition replaced by the relaxed stability eigenvalue 
condition. Then, for sufficiently small e, the problem is a contraction, in a suitable H- 
norm, for the nontrivial part v of the solution of (1.1) and — > const, when t — > oo. 
Thus, the system (1.1) is stable in this generalized sense. 
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